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SUMMARY 

In  design  problems  for  fusion  reactor  structures  and  other 
magnetically  loaded  devices,  the  field  and  the  magnetic  forces  are 
usually  calculated  for  static  conditions  only.  The  structural  analysis 
is  then  handled  as  a  separate  problem.  Recently,  more  attention  is 
being  given  to  the  transient  field  problems,  but  little  work  has  been 
done  in  understanding  the  nature  of  the  dynamic  magnetic  forces  and  the 
mutual  interactions  between  the  field  and  the  structure.  This  paper 
presents  a  study  of  the  magnetic  forces,  the  mutual  interaction  between 
fields,  and  the  induced  vibrations  of  long  conducting  plates.  The 
finite  element  method  is  used  for  analysis.  Some  numerical  results  are 
presented. 

INTRODUCTION 

Magnetic  forces  result  from  the  interactions  between  the  external 
field  and  the  induced  current  in  a  conducting  body,  assumed  here  to  be 
nonferromagnetic.  The  phenomenon  of  currents  induced  by  time-varying 
magnetic  fields  is  governed  by  Faraday's  induction  law,  and  is  called  the 
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"eddy  current  problem1'  in  electromagnetism.  The  problem  is  a  three- 
dimensional  one  for  conducting  plates  with  finite  thickness.  When  the 
thickness  of  the  plate  is  small  compared  to  the  penetration  depth  of 
the  magnetic  field,  the  induced  current  density  across  the  thickness  of 
the  plate  is  approximately  uniform.  A  stream  function  representation  of 
the  current  may  then  be  used  to  reduce  the  problem  to  a  two-dimensional 
one.  It  is  usually  assumed  that  the  externally  applied  magnetic  field 
is  not  affected  by  the  induced  current.  The  total  magnetic  field  may  be 
decomposed  into  an  externally  applied  part  and  a  self-induced  part  due 
to  the  generation  of  the  eddy  current.  The  self-induced  part  may  be 
determined  by  the  Biot-Savart  law,  which  involves  integration  of  field 
quantities  over  the  entire  volume  of  the  plate.  An  integro-differential 
equation  may  then  result  for  the  determination  of  the  stream  function. 

The  integral  terms  represent  the  flux  linkage  of  the  eddy  current  density 
over  different  part  of  the  plate  and  are  called  the  "non-local"  terms 
below.  For  infinitely  long  plates,  a  one-dimensional  integro-differential 
equation  may  be  derived  for  the  calculation  of  the  induced  eddy  current. 
Results  of  a  study  for  steady-state,  harmonic  currents  in  long,  rigid 
plates  have  been  presented  in  Reference  1. 

EQUATIONS  AND  FINITE  ELEMENT  FORMULATION 

The  basic  equations  for  the  coupled  theory  for  linear  nonferro¬ 
magnetic  plates  are  the  Maxwell's  equations 

7  x  H  =  J  7  •  B  *  0 

7xE  =  -^rB  7  •  E  =  0 

supplemented  by  Ohm's  law 

J=aLE  +  ^-uxB] 


(2) 
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and  the  linear  equation  of  motion  for  the  conducting  plate 

2 

DV4w  +  ph  — y  =  F  +  n  •  V  x  C 

n 


where 


F  = 


C  = 


h/2 

n  •  (J  x  B)  dz 
-h/2  ~ 

h/2 

n  x  (J  x  B)  z  dz 
-h/2  " 


(3) 


(4) 


in  which  E,  H,  B,  and  J  are  the  electric  field  intensity,  the  magnetic 
field  intensity,  the  induction,  and  the  current  density,  respectively, 
u  is  the  displacement  field  and  w  the  transverse  deflection  of  the  mid¬ 
surface  of  the  plate,  n  is  the  unit  normal  vector  to  the  mid-surface. 
a  is  the  electrical  conductivity,  p  the  mass  density,  h  the  thickness  of 
the  plate,  and  D  the  bending  rigidity  of  the  plate.  (A  complete  list 
of  symbols  used  is  given  in  the  Appendix). 

Using  the  stream  function  representation  for  the  current,  I  = 

7x  (<j>n)  in  which  I  =  hJ,  and  following  a  procedure  suggested  by  Moon  [2] 
for  the  calculation  of  the  self-field,  one  obtains  the  following  equations 
for  long  conducting  plates  (Fig.  1): 


djjj  3  r  f 

T7  -  K  CoV  -  IT  j 


=  ah 


3B 

r  z  .  _L  (lit 
L3t  3x  '3t  VJ 


*(£)  -I 

(S-x)2  +  lh2 


(5) 
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+  D 


34w 
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in  which  the  superscript  0  represents  the  applied  (external)  portion  of 
the  total  field.  The  time-dependent  terms  on  the  left  hand  side  of  Eq. 
(5)  represent  the  nonlocal  effect  of  the  self-field.  In  both  equations 
the  coupling  effects  appear  on  the  right-hand  side  only. 

The  finite  element  Galerkin  method  is  used  to  solve  the  coupled 
equations  (5)  and  (6).  \ p  is  approximated  globally  and  locally  by 

piece-wise  linear  models 

G  2  r 

ip  =  E  M.  vp.  and  <jj  =  Z  N.  tk  (7) 

k=l  K  K  k=l  K  K 

in  which  G  is  the  total  number  of  nodal  points,  Mk  are  the  global 
interpolation  functions  generated  from  the  local  linear  element  shape 
functions  with  the  superscript  E  denoting  the  Eth^ element,  w  is 
approximated  by  the  usual  cubic  model. 

6  E 

w  =  Z  w  (8) 

k=l  k  k 

The  following  set  of  linear  algebraic  equations  for  each  element  then 
results  from  Eqs.  (5)  and  (6): 


"  Sjk  \  Pjk  +  Qjk  *k  =  Rj  (9) 


in  which 


(  dNj  dN^ 

Sjk  =  I  "dJ  "d7  dx 


Pjk  =  <*o  f  Nk  dx 
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c  om  E 

4  ■  7T  J  ¥«>  “j<5)  « 

0 

Rj  '  f  "5  CB°  -  £  dx  (10) 

where  the  weighting  function  is 


t  r 

Both  W.(5)  and  Q..  are  integrated  analytically.  The  integra- 
J  J  K 

tions  for  are  carried  out  over  the  entire  plate  and  result  in  a  full, 
usually  unsymmetric,  matrix  [3]. 

In  matrix  notation,  the  equations  of  motion  are 

[M]{w)  +  [K]{w}  =  {F}  (12) 

in  which 

{F>  =  -  f  (C)  ||  B°  dx  (13) 

E 

The  global  matrix  form  of  Eq.  (9)  is 

[A]{*}  -  [S]W  -  TR>  H4) 

in  which 

[A]  -  -  [P]  +  [Q]  (15) 

In  the  general  case,  [A]  is  full  and  nonsymmetric. 

The  transient  equations  (14)  and  (12)  are  integrated  using  the 
following  scheme. 


*♦  . 
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<3t  W  -  d-e)[s]}  (*}tt4t 

a  ^  +  9CS]}  CY}t  +  9(R)t  +  (1-9)  {R;t+4t  06) 

and 

m  *  w>  («t+4t 

2 

=  {F}t+At  -  [M]  {{w}t  +  At{w)t  +  {«}t}  (17) 

Equations  (16)  and  (17)  form  a  weakly  coupled  set  of  equations.  The 

3  *0 

forcing  function  {R}  in  (16)  contains  terms  proportional  to  ^(wBx). 

Its  magnitude  is  small  compared  to  the  total  induced  eddy  currents.  The 
computation  is  performed  by  two  more  or  less  independent  routines.  The 
solution  to  the  coupled  system  is  advanced  by  sequentially  executing 
these  two  routines.  The  coupling  term  in  Wt+At  in  Eq.  (16)  is  esti¬ 
mated  by  temporal  extrapolation. 

EXAMPLES 

Numerical  simulation  has  been  conducted  for  the  magnetically 
induced  vibrations  of  thin  aluminum  plates  (conductivity  a  =  3.8  x  107 
1/ohm-m  and  wave  speed  cy  =  /e7p  =  5090  m/sec).  Because  of  the  widely 
different  time  response  characteristics  of  the  electromagnetic  and 
mechanical  fields,  two  levels  of  investigations  have  been  carried  out. 

To  study  the  effects  of  the  magnetic  coupling  of  the  eddy  current  density 
on  the  distribution  and  time  variation  of  current  and  force,  the  plate 
is  first  assumed  to  be  rigid.  For  the  arrangement  shown  in  Fig.  1(c), 
with  l  =  485  mm,  a  =  94.31mm,  b  =  314  mm,  e  =  9.922  mm  and  h  =  2.381  mm, 
the  transient  analysis  for  a  half-sine  driving  pulse  in  a  double-wire 
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exciting  coil  gives  the  results  shown  in  Fig.  2.  The  effect  of  the 
nonlocal  terms  is  apparent  and  includes  both  a  phase  shift  and  a  nonzero 
average  pushing  force,  as  predicted  in  the  study  for  steady-state, 
harmonic  current  cases  [1], 

The  coupled  magnetomechanical  effect  is  considered  for  the 
arrangement  shown  in  Fig.  1(b),  with  i  =  427  mm,  a  =  334.6  mm,  h  =  2.11  mm 
and  e  =  4h  =  8.44  mm.  The  exciting  wire  is  purposely  placed  above  the 
nodal  point  of  the  second  vibration  mode  of  the  cantilever  plate. 

Because  of  the  much  sharper  variation  of  the  electromagnetic  field 
variables,  more  field  elements  than  beam  elements  are  used  in  the  analysis 
The  meshes  of  the  two  subsystems  are  made  conformable  by  dividing  each 
beam  element  into  several  equal-length  field  elements.  The  coupling 
terms  in  the  two  subsystems  are  evaluated  by  interpolation  and  numerical 
integration. 

To  have  meaningful  coupling  between  the  electromagnetic  and 
mechanical  subsystems,  a  3ms  pulse  duration  in  the  exciting  wire  is 
chosen,  half  the  period  of  the  third  vibration  mode  of  the  beam-plate 
(T^  =  6ms,  At  =  0.075ms).  The  results  of  the  analysis  are  shown  in 
Figs.  3  and  4  for  a  driving  current  of  1°  =  500  Amp.  The  u  x  B  term 
("two-way  coupling")  demonstrates  itself  as  a  damping  effect  and  has 
its  maximum  influence  in  the  region  close  to  the  exciting  wire.  The 
time  histories  of  eddy  current,  force,  and  transverse  deflection  at  a 
point  close  to  the  exciting  coil  are  shown  in  Figure  4.  The  force  plot 
chows  that  a  pulling  force  can  be  generated,  and  the  displacement  plot 
indicates  that  the  energy  transfer  into  the  mechanical  subsystem  is 
reduced  because  of  the  u  x  B  effect.  The  nonlocal  effect  is  significant 
and  generates  a  net  pushing  force.  The  third  mode  component  also  shows 


Deflection  (mm)  Magnetic  force  (N/m)  Eddy  current  density  (K  amp/ 


Figure  3.  Nonlocal  predictions  of  eddy  current  density,  mag¬ 
netic  force,  and  transverse  displacement  at  t  = 

0.9  ms  for  cantilever  plate  with  single-wire 
exciting  coil . 


Transverse  displacement  (nun)  Magnetic  force  (N/m)  Eddy  current  density  (Kamp/m) 
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Figure  4.  Nonlocal  predictions  of  eddy  current  density, 
magnetic  force,  and  transverse  displacement 
at  x  s  33  cm  for  cantilever  plate  with  single 
wire  exciting  coil . 
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up  clearly  in  the  displacement  plot,  superimposed  on  the  first  mode 
response  (T^  =  105ms). 

CONCLUSION 

A  coupled  electro-magneto-mechanical  problem  has  been  studied  for 
the  magnetically  induced  vibration  of  thin  nonferromagnetic  conducting 
plates.  Some  transient  analysis  results  have  been  presented  for  the 
pulsed  current  problem.  The  nonlocal  and  u  x  B  couplings  have  been 
investigated.  The  results  show  that  the  nonlocal  effect  is  essential 
for  the  calculation  of  force,  and  that  the  u  x  B  term  acts  as  a  damping 
effect  in  the  coupled  problem.  Experience  so  far  has  been  limited  to 
simple  pulsed  currents.  Work  in  progress  includes  experimental  verifi¬ 
cation,  further  numerical  studies,  and  extension  to  geometrically  nonlinear 
effects  of  deformation.  Moreover,  cyclically  pulsed  exciting  currents 
need  to  be  examined  as  a  possible  source  of  structural  instability. 
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APPENDIX:  LIST  OF  SYMBOLS 
[A]  electromagnetic  field  matrix 
B  induction 

D  bending  rigidity 

E  electric  field  intensity 

F  magnetic  force 

h  thickness  of  plate 

H  magnetic  field  intensity 

1°  driving  current  in  wire 

I  current  per  unit  length 

J  current  density 

[K]  stiffness  matrix 

l  width  of  plate 

M^  global  shape  functions 

[M]  mass  matrix 

n  unit  normal  vector  to  the  mid-surface 

N^,C^  linear  and  cubic  element  shape  functions,  respectively 
[P]»[Q].[S]  electromagnetic  field  matrices 
u  displacement 

w  transverse  deflection  of  plate 

E 

W.  weighting  functions 

J 

x,z,C  orthogonal  cartesian  coordinates 

u3  permeability  of  free  space 

p  mass  density 

a  electric  conductivity 

p  stream  function  for  eddy  current 

0  superscript  0  denotes  externally  applied  quantities 
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